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We analyze the influence of dipole-dipole interactions in an electromagnetically induced trans- 
parency setup for a density at the onset of cooperative effects. To this end, we include mean-field 
models for the influence of local field corrections and radiation trapping into our calculation. We 
show both analytically and numerically that the polarization contribution to the local field strongly 
modulates the phase of a weak pulse. We give an intuitive explanation for this local field induced 
phase modulatton and demonstrate that it distinctively differs from the nonlinear self-phase modu- 
lation a strong pulse experiences in a Kerr medium. 
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I. INTRODUCTION 



Electromagnetically induced transparency (EIT) 
stands out as one of the most useful coherence and 
interference phenomena (see P, and references 
therein). Current research focuses on dilute samples 
with iVA'^ <C 1 {N density, A transition wavelength), 
in which the atoms essentially act independently. Ex- 
perimentally, an important reason for the restriction 
to low densities is detrimental decoherence induced, 
e.g., by atom collisions. This density restriction applies 
in particular to hot atomic vapors. More promising 
in this respect are ultracold gases [sl-fioj. in which 
atomic collisions are much less frequent, leading to 
greatly improved coherence properties. However, most 
experiments still operate in the regime of a dilute gas, 
where cooperative effects do not play a role. But with 
recent advances in preparation techniques, now trapped 
ultracold gases at densities up to 10^'' — 10^^ atorns/cm'^ 
with a linear extend of the gas in the fim range |17H19| 
have been reported. With such setups, dense gas light 
propagation experiments at the onset of cooperative 
effects seem within reach. 

From the theoretical side, it is well-known that in the 
high density regime (A^A^ 3> 1), a rigorous treatment 
of cooperative effects in general gives rise to an infinite 
hierarchy of polarization correlations [20l - [2^ . Calcula- 
tions of this type are demanding, such that concrete re- 
sults for more complex model systems such as EIT have 
not been found yet. However, at the onset of coopera- 
tive effects, the infinite hierarchy can be truncated. A 
study of optical properties of a dense cold two-level gas 
found that in leading order of the density, three correc- 
tions with respect to a dilute gas occur [23]. First, the 
microscopic field Emic driving the individual atom is no 
longer the externally applied field E^xt alone, but rather 
is corrected by the mean polarization P of the neighbor- 
ing atoms. This local field correction (LFC) is described 



by the well-known Lorentz-Lorcnz formula [24 
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The second correction in leading order of the density is 
due to the quantum statistics of atoms, and becomes rel- 
evant close to the phase transition to a Bose-Einstein 
condensate (BEC). In this regime the atomic de Broglie 
wavelength is of the order of the electronic transition 
wavelength such that the medium's refractive and ab- 
sorptive properties are enhanced by quantum mechanical 
exchange effects. In contrast, sufficiently above or below 
the transition point, quantum statistical corrections are 
small [23]. The third correction arises due to the leading 
order of multiple scattering, which can be interpreted as 
a reabsorption of spontaneously emitted photons before 
they leave the sample. In general, the reabsorption gives 
rise to a complicated dynamics because of the many path- 
ways a photon can take between the different atoms in 
the medium. But to a first approximation, this radiation 
trapping can be modeled in a mean field like approach 
by suitably chosen additional incoherent pump rates be- 
tween the ground states and the excited state [25l - [28t . 

These results suggest that at the onset of cooperative 
effects, and at parameters away from the phase transi- 
tion to BEC, the dominant cooperative corrections arise 
on a mean-field level, such that a macroscopic treatment 
is meaningful. Similar results were obtained in a re- 
cent calculation which for a two-level system explicitly 
compares a microscopic treatment based on a multiple- 
scattering approach including n-particle correlations with 
a macroscopic treatment based on LFC [2^ . It was found 
that with densities of few particles per wavelength cubed 
(A^A'^ w 5), the two methods agree well, while at higher 
density (A^A'^ « 124) deviations occur, which however 
due to the numerical complexity could only be stud- 
ied for rather small sample sizes. So far, LFC alone 
has been studied in a number of systems (soj . includ- 
ing self-induced transparency for short and ultrashort 
pulses dll, [11] , coherent population trapping [sl] and 
lasing without inversion [34j. More recently, LFC was 
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FIG. 1: (Color online) EIT level scheme with probe field Q31 
(red dashed arrow) and control field ^32 (blue solid arrow). 
The green wiggly arrows indicate the radiative decay 731 , 732 
and the incoherent pump rates R31, R32, while 7deph, 7s de- 
note the ground state dephasing and the ground state popu- 
lation transfer. 



studied in the context of atomic gases with a negative 
refractive index (NRI) [35l439l] . LFC effects have also 
been exploited to spectroscopically resolve the hyperfine 
structure of the Rb D2 line in a dense gas [i^ . 

But surprisingly, cooperative effects in EIT have re- 
ceived only little attention. In [25l,[2^, effects of radiation 
trapping on EIT were studied. The propagation of two 
non-adiabatic propagating pulses is considered in (4ll |. 
Some effects such as a modification of the group velocity 
and a phase modulation are reported, but only numerical 
results are given without clear interpretation [Z^l . 

The recent theoretical proposals for NRI requiring high 
density already motivate a better study of the regime of 
higher densities in atomic gases. But also several as- 
pects of EIT itself could benefit strongly from such an 
analysis. For example, the group velocity reduction in 
EIT and the related spatial compression of the pulse are 
directly proportional to the density Also, the time 
delay bandwidth product, a figure of merit for the over- 
all performance of a slow light system, depends on den- 
sity [13 . Finally, a higher density could facilitate a more 
efficient coupling of light and matter, e.g., for applica- 
tions in quantum information science p^ j. 

Therefore, in this paper, we study light propagation 
in an EIT medium at the onset of cooperative effects, 
and reveal and interpret the underlying physical mecha- 
nisms. We show that LFC leads to a phase modulation 
of a probe pulse at densities achieved in current exper- 
iments. This phase modulation is distinctively different 
from nonlinear self-phase modulation, as it leads to a 
linear frequency chirp across the whole probe pulse, and 
does not require high probe intensity. Our main aim is 
a physical understanding of the cooperative effects. For 
this, we include LFC, but assume parameters away from 
the crossover point to a Bose-Einstein condensate, such 
that corrections due to quantum statistics can be ne- 
glected. Similarly, our semi-classical treatment of the 
light fields does not allow to reveal non-classical proper- 



ties of the scattered light [i^. We model the onset of 
multiple scattering via an incoherent pump rate propor- 
tional to the excited state population. Furthermore, we 
take into account a ground state decoherence due to elas- 
tic collisions as well as a ground state population transfer 
due to spin-exchange collisions. In addition to a full nu- 
merical solution of the equations of motion (EOM), we 
derive analytic solutions for the propagation dynamics 
of a slow light pulse and the medium polarization, which 
enable us to provide an intuitive explanation of the phase 
modulation in terms of the energy exchange with neigh- 
boring atoms. 

The paper is organized as follows. In the next section 
(Sec. |TT| the theoretical model is described. We present 
the full set of equations of motion (EOM) (Sec. Ill Al) . 
discuss our framework for taking into account radia- 
tion trapping (Sec. IIIBI) . and introduce LFC into the 
EOM (Sec. Ill C|) . In Sec. IIIII we present our results. 
First, we give an analytic expression for the susceptibil- 
ity (Sec. nil A[) . Then, we use this expression to derive 
an analytic solution for a propagating probe pulse in- 
cluding the LFC-induced phase modulation (Sec. IIIIBI) . 
In Sec. nil CI we compare this analytic solution to the 
numeric solution of the full set of equations of motion 
and calculate the expected probe field absorption due to 
radiation trapping and ground state population trans- 
fer for a range of different possible experimental param- 
eters. In Sec. IIIIDI we explain the physical origin of 
the phase modulation in terms of the energy exchange of 
neighboring atomic dipoles and contrast it with the self- 
phase modulation a strong pulse acquires in a nonlinear 
medium. Finally, in Sec. IIVI we draw some conclusions. 



II. THEORETICAL MODEL 

A. Maxwell-Schrodinger equations 

We start from the Hamiltonian for an EIT level 
scheme 0] as shown in Fig. [T] In a suitable interaction 
picture it is given by 

H = ^ hiAA22 + A31A33) (2) 
- ^ (f^aiAai + f^32^32 +h.c.) , 

where A = A31 — A32 is the two-photon detuning, 
Aij = \i){j\ are the atomic operators, and ^32, 
A31, A32 are the Rabi frequency and detuning of the 
probe and control field, respectively. Note that without 
application of the LFC, the Rabi frequencies contain mi- 
croscopic fields at the position of the respective atoms. 
After performing LFC, in our numerical calculations, the 
microscopic electric fields are replaced by the externally 
applied fields as discussed in Sec. Ill CI 

In order to account for the various incoherent pro- 
cesses, we choose a master equation approach |4(|], and 
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the EOM for the atomic density operator p reads 
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(3a) 



^^([pA3„A,3]+h.C.) 



^^([A3„[A,3,p]]+h.C.) 

3 = 1 

-^([A21,[Ai2,p]]+h.C.) 
- Ideph {A22pAii + h.C.) . 

Here, 73^ [j G {1,2}) denote the radiative decay rates, 
and and R^j are incoherent pumping rates between the 
ground and the excited state. 7^ is the ground state 
population transfer rate, and ^deph is the ground state 
dephasing. 

The dynamics of the probe and control fields are each 
governed by a wave equation. In slowly varyirig_envelope 
approximation (SVEA), these have the form [43] 



1 



dt 



ik 

£33 = — P3j 

So 



j€{l,2}, (3b) 



where c is the speed of light in vacuum, S^j is the envelope 
of the probe and control field, and P^j is the probe and 
control field polarization. 

The combination of Eq. pap and Eqs. pb|) form the full 
set of Maxwell-Schrodinger EOM that govern the spatio- 
temporal dynamics of the system. In our model, its solu- 
tion becomes more challenging because of the additional 
nonlinear terms introduced due to radiation trapping and 
LEG. We will discuss both aspects separately and in more 
detail in the following. 



dynamics. In an EIT system, the excited state popula- 
tion in principle should remain small. Especially under 
adiabatic conditions when the system is in the dark state 
at all times, the destructive quantum interference under- 
lying EIT strongly suppresses excitation events. How- 
ever, in experiments a dephasing of the ground state co- 
herence as well as spin exchange due to inhomogeneous 
magnetic fields and atomic collisions can lead to an effec- 
tive pump rate into the bright state and subsequent exci- 
tation. This type of dynamics is modeled by the different 
incoherent pump rates in the master equation Eq. (j3al) . 
In Sec, nil Cl we use the numerical solution of the full set 
of EOM including this dynamics to calculate the probe 
field absorption for a range of different values for r^/re 
as well as 7s. 



Local field corrections 



In the master equation (|3a|) the Rabi frequencies de- 
pend on the microscopic probe and control field, such 
that we have to replace them by their macroscopic coun- 
terparts using Eq. ([1]). Since the mean medium polar- 
ization P can be expressed in terms of a density matrix 
element, this leads to nonlinear EOM. Expanding the 
nonlinear EOM up to linear order in the external probe 
field leaves us with just two EOM, one for the probe field 
coherence psi and one for the Raman coherence p2i, 

XX X 
dtPsl = - r3iP3i + -1^31 + -£731^31 -I- -f^32P21 , (5a) 



dtP21 = - r2i/921 + 2^32/031 



(5b) 



From now on, after performing the LEG, il^i and r232 
denote the respective Rabi frequencies of the external 
fields, and we have defined 



B. Radiation trapping 

A higher density of the atoms leads to reabsorption 
and multiple scattering of spontaneously emitted pho- 
tons. We model this by incoherent pump fields R31 and 
i?32 with an intensity proportional to the population of 
the excited state. The number of incoherent photons in- 
side the sample is also infiuenced by its geometry. Along 
the lines of 27[, this is taken into account by a pumping 
rate ra due to atomic decay and a photon escape rate r^. 
Altogether, R31 and i?32 are then given by 



^31 = 731 — P33 

D _ ''a 
^32 — 732 — P33 



(4a) 
(4b) 



Obviously, the ratio ra/r^ should also depends on the 
atomic density. For high density it will eventually ap- 
proach unity such that even a small fraction of excited 
state population can lead to a significant change of the 



r3i =^ - «A3i, 

r21 =7dec - iA, 

7 =731 +732 +7s, 

Idee =ldeph + 7s • 



(6a) 

(6b) 
(6c) 
(6d) 



Due to LEG, a new term arises in Eq. (I5al) as compared 
to the low-density case which is proportional to the di- 
mensionless parameter 



L 



4^ 



(7) 



It is a measure for the strength of LEG, and the prefactors 
ensure that L = 1 corresponds to a density where the 
LEG induced frequency shift in a two-level atom is equal 
to half the natural linewidth. Formally, the new term 
can be interpreted as a frequency shift in an EIT system 
as well and can be included into the probe field detuning 
A31 = A31 -I- L731/2. However, this frequency shift does 
not infiuence the two-photon detuning A. 
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FIG. 2: (Color online) Real part (blue dashed line) and imag- 
inary part (red solid line) of the EIT susceptibility given by 
Eq.[8]for (a) negligible LFC (L = 10"^) and (b) strong LFC 
(L = 4). 



III. RESULTS 

A. Susceptibility with LFC 

The EOM ([S]) are already linear in the probe field and 
solving for the steady state, we can derive the suscepti- 
bility X, 



3L^(A + ndec) 



A31A 



- i(A3i7dec + A^) 



(8) 



With an analytic expression at hand we can easily pin- 
point the effect of LFC. We find that a reshaping of the 
EIT transparency window takes place [s^. To illustrate 
this, in Fig. [2] we show the real and imaginary part of x 
for two different sets of parameters. In Fig. [H^a) we set 
L = 10~^ which results in negligible LFC and the stan- 
dard form of the EIT susceptibility is recovered. With 
L = 4 in Fig. [2][b) a strong reshaping of the transparency 
window due to LFC is found. In both cases the control 
field Rabi frequency is ^-yi — 27 and the control field 
detuning is A32 — 0. 



B. Analytic solution for a light pulse 

To analyze the infiuence of LFC on the propagation 
dynamics of a light pulse, we expand Eq. ([8]) around the 
center of the transparency window. From 



(9) 



we find the frequency dependent wave number h(y)) up 
to second order in the probe field detuning. The solution 



for the positively rotating component of the probe field 
in Fourier space then follows from 

(z, cj) = (0, w) exp[ifc(w)z] , (10) 
where £:(+)(0,a;) is given by the initial condition and 



fc(w) « fco + i 



!V]^ + ^ + fco(^^,+/3,)A2^. (11) 



We neglected terms suppressed by a factor of ^^^ecj^i, 
since ^^dec ^ ^^§2 is required for low absorption, and of 
higher order in A31. Each term in Eq. (Ilip can be clearly 
interpreted, fco is the wave number of the undisturbed 
carrier wave. The second teriu describes the decay due 
to the ground state decoherence ^dec where 



3X73 ICJQ 

^32 



(12) 



is the group index and is the probe field transition 
frequency. The third teriu leads to the reduced group 
velocity Vg = c/(l + rig). The fourth term is quadratic 
in A31 and thus associated with a change of width in the 
Fourier transformation of a Gaussian. The imaginary 
part proportional to 



/3i 



6Jv773i 

^32 



(13) 



leads to a broadening of the temporal width due to the 
finite spectral width of the transparency window. Finally, 
the real part proportional to 



/32 = 



^32 



(14) 



results from LFC. Formally, it leads to an imaginary part 
in the temporal width which corresponds to a phase mod- 
ulation of the pulse. When we assume a Gaussian pulse 
shape for the initial condition, the solution in the time 
domain can be obtained by Fourier transformation. Con- 
sidering only the envelope defined by 

E{z, t) — ^^{z, t) exp[i(fcoz — ujot)] + c.c. , (15) 



we find 



£iz,t) 



OQ — exp 
a 



-Ide 



Z 

'v„ 



2d 



S-2 



(16) 



£q and a are the initial amplitude and temporal width. 
The LFC luodified width after propagating a distance z 
is 



2kozil3i - 1^2) 



(17) 



The phase modulation cair be well approxiiuated by the 
parabola 



l32koz 



1-it-z/vgf /a' 



(18) 
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FIG. 3; (Color online) Probe field absorption in dependence 
of the radiation trapping rate Taj re and the ground state 
transfer rate 7s. The shading indicates the peak value of the 
propagated probe field amplitude ^31(2) relative to the initial 
amplitude ^31(0). The three contour lines indicate parame- 
ters at which 50%, 25% or 1% of the initial amplitude leave 
the sample. The results are obtained using the numerical so- 
lution of the full set of EOM. All other parameters remain 
fixed and are chosen as discussed in Sec. HIT CI 

C. Numeric solution 

To compare our analytic solution to results obtained 
numerically, we chose parameters consistent with recent 
ultra cold gas setups |17l - [l9{ . We assume a medium with 
density iV w lO^^cm"^ and a length of z « 40/xm. This 
corresponds to N)? w 50 (A = 795nm), and with a con- 
trol field f232 = 27, a probe pulse with initial width of 
a = 20/7 propagates a distance of z 150vg/7. Besides 
some broadening due to the finite spectral width of the 
transparency window, the probe pulse is attenuated di- 
rectly by the ground state decoherence ^dec (see Eq. . 
It also suffers absorption because of the reabsorption of 
incoherent photons. As discussed in Sec. IIIBl this signif- 
icantly depends on the ground state population transfer 
7s and the geometry of the sample represented by /re . 

Prom results of a recent experiment with Rubidium 
atoms of density iV « 2 x lO^^cm"'^ in an anisotropic 
trap we estimate ^dec to about 100 kHz [T7l - [l9| . From 
radiation trapping experiments in hot and ultracold 
gases [13, , we further estimate that the ratio of the 
photon pump and escape rate is about rajr,, = 0.99. The 
rate of spin exchange collisions depends on the relative 
ground state energies and the fact that spin exchange 
collision become more likely if the overall angular mo- 
mentum is conserved. From (48j we find that a typical 
value in ultracold gases is given by 

7s = X lO^^cm^s"^ , (19) 

such that for a density of A^ ?» lO^^cm"^ we estimate at 
a value of about 7^ = 1 kHz. To account for different 
experimental situations, and to obtain insight in the ef- 
fect of radiation trapping and ground state decoherence, 
we use the full set of EOM to numerically propagate a 



probe pulse and to calculate the probe field absorption 
for a range of values for 7s as well as Va/re- The result 
is shown in Fig. |3l The plot shows the amplitude of the 
propagated pulse relative to the initial pulse. The three 
contour lines indicate parameters for which a fraction of 
50%, 25% or 1% of the initial light field amplitude leaves 
the sample, respectively. For the estimated values of 7dec, 
7s, and Vajr^ which we assume for our further calcula- 
tions, the probe field f^ai is reduced to about 50% of its 
initial value. 

Fig. |4] analyzes the pulse after passing through the 
medium. In subfigure|31[a), we show the propagated pulse 
together with the LEG induced phase modulation 0lfc (0 
and the corresponding instantaneous frequency defined 
by w(i) = cjQ — dt(f>i,Fc{t)- The numerical solution is vir- 
tually indistinguishable from the analytical one, and we 
only show the analytical solution. We see that the phase 
modulation has a negative parabola-like shape such that 
the instantaneous frequency is approximately linear over 
the total extend of the pulse. 

D. Physical origin of the phase modulation 

To explain the physical origin of the phase modulation, 
we explicitly calculate the relevant parts of the polariza- 
tion using the relation 

P^+Hz,uj) = eox{uj)E^+\z,uj) . (20) 

Considering only the real part of x up to quadratic order 
in A31 and Fourier transforming it into the time domain, 
we can distinguish two contributions, 

Pl^+\z,t) = ^ [idt£{z,t)]exp[i{koz - ujot)] , (21a) 

P[t^{z, t) - eo/?2 [i^d^Siz, t)] exp[z(A:oz - ujot)] . (21b) 

The first contribution stems from the part linear in A31 
and leads to the change of group velocity. The second 
contribution is due to the part quadratic in A31 which 
is related to the LEG induced phase modulation. In 
Fig. mjb), we show the pulse together with these two 
contributions. In the first half of the pulse, E{t) is ahead 
in phase by 7r/2 compared to Poit), which indicates that 
energy is transferred from the pulse to the polarization 
Po{t). In the second half, the pulse is delayed by tt/2, and 
energy is transferred back from the polarization Po(t) to 
the pulse. This energy exchange effectively reduces the 
group velocity of the pulse. Similarly, we can understand 
how the interaction of atoms with the collective dipole 
field of their neighbors proportional to Po(t) induces an 
additional polarization PLFc(i)- Before t = —a, the po- 
larization component Po(t) is 7r/2 ahead in phase com- 
pared to PLFc(i), whereas at —a < i < 0, it is delayed 
by 7r/2, again leading to an energy exchange. The same 
exchange takes place again for < < < cr and t > a. 
While the additional polarization Plfc(0 is induced ex- 
actly by the same mechanism as Po{t), its back action on 
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FIG. 4: (Color online) (a) Real part (blue dashed line) and 
imaginary part (red solid line) of a Gaussian pulse with LFC 
induced phase modulation. Both parts are scaled to arbitrary 
but equal units. The parabola (black dash-dotted line) shows 
the time dependent phase shift in radian measure (right axis) 
and the dashed straight line is the corresponding instanta- 
neous frequency. It exhibits a linear frequency chirp with 
slope Olfc over the total extend of the pulse, (b) Envelope 
and carrier wave of the pulse (black upper line), the polar- 
ization induced directly by the pulse (blue middle line), and 
the additional polarization induced by dipole-dipole interac- 
tions (red lower line). The amplitudes are not drawn to scale 
and the carrier wave length has been strongly exaggerated 
to make phase relations clearly visible, (c) The nonlinear 
self-phase modulation (black dash-dotted line) of a Gaussian 
pulse. The frequency chirp is approximately linear only in the 
pulse center (black dashed straight line) with slope qnsm. 



the probe pulse is different. At t < —a and t > cr, PlpcI*) 
and E{t) have opposite phase. This opposite phase has a 
dragging effect on E{t), and reduces its phase. In the cen- 
tral part of the pulse {—a < t < cr), Plfc(^) is in phase 



probe pulse is different. Ait < —a and t > cr, Plpc 
I E(t) have opposite phase. This opposite phase ha; 



with E{t). This has a pushing effect on E{t), and in- 
creases the phase of the pulse. This interpretation agrees 
with the phase modulation obtained from the calculation 
as shown in Fig. IH^a). We thus conclude that energy is 
exchanged between the atomic dipoles and the field of 
neighboring dipoles in exactly the same way as between 
the atomic dipoles and the external field E{t). But the 
two polarization components act differently on the probe 
pulse, leading to the group velocity change and the phase 
modulation, respectively. 

Finally, we compare the LFC induced phase modu- 
lation with nonlinear self-phase modulation (NSM) in 
a medium with an intensity dependent refraction. The 
NSM modulation is 



0NSM(i) = n2l{t)koz, 



(22) 



where 712 is the intensity dependent index of refraction 
and I{t) the intensity profile of the pulse. In Fig. \M,c) 
we show 0NSM (t) together with the corresponding instan- 
taneous frequency for a Gaussian pulse. We see that the 
front of the pulse experiences a red shift whereas the back 
experiences a blue shift with an approximately linear fre- 
quency chirp in the center. Comparing the two chirps, 



2n2/ofcoz 
2/32 fco^ 



(23) 
(24) 



we find that in the LFC case, 71,2/0 is replaced by /32/o'^- 
Thus, the LFC modulation does not require a large inten- 
sity, and is approximately linear over the total extend of 
the pulse since it depends on the strength of the dipole- 
dipole interaction. This strength is given by in an 
EIT system and can be infiuenced by the density and 
the control field strength ri32. 



IV. CONCLUSION 

In summary, we studied light propagation in an EIT 
system for densities at the onset of cooperative effects. 
For this, we amended the standard Maxwell- Schrdinger 
equations for a three-level A-type EIT setup on a mean- 
field level by local field corrections of the applied fields, 
and by a model for radiation trapping based on inco- 
herent pump rates between the ground and the excited 
states. As main result, we found that the phase of a 
Gaussian probe pulse is strongly modulated by the polar- 
ization contribution to the local field. We gave a physical 
interpretation for the underlying mechanism leading to 
the phase modulation and showed that it is distinctively 
different from the case of non-linear self phase modula- 
tion, as the resulting frequency chirp is linear across the 
whole pulse, and does not depend on the pulse intensity. 
Based on the numerical solution, we further studied the 
effects of decoherence and radiation trapping. Within 
our model, this allows to estimate upper limits for these 
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incoherent effects such that the pulse survives a sufficient observable, 
propagation length for the phase modulation to become 
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